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In this talk we present data analysis methods for reconstructing the mass and 
couplings of Weakly Interacting Massive Particles (WIMPs) by using directly 
future experimental data (i.e., measured recoil energies) from direct Dark Mat- 
ter detection. These methods are independent of the model of Galactic halo 
as well as of WIMPs. The basic ideas of these methods and the feasibility and 
uncertainties of applying them to direct detection experiments with the next 
generation detectors will be discussed. 
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1. Introduction 

Weakly Interacting Massive Particles (WIMPs) x arising in several ex- 
tensions of the Standard Model of electroweak interactions with masses 
roughly between 10 GeV and a few TeV are one of the leading candidates 
for Dark Matter ■ . Currently, the most promising method to detect differ- 
ent WIMP candidates is the direct detection of the recoil energy deposited 
by elastic scattering of ambient WIMPs on the target nuclei ' . The dif- 
ferential event rate for elastic WIMP-nucleus scattering is given by-*^: 



dR I poco 



fiiv) 



dv . (1) 
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Here R is the event rate, i.e., the number of events per unit time and unit 
mass of detector material, Q is the energy deposited in the detector, po is 
the WIMP density near the Earth, ctq is the total cross section ignoring 
the form factor suppression, F{Q) is the elastic nuclear form factor, fi{v) 
is the one-dimensional velocity distribution function of incident WIMPs, v 
is the absolute value of the WIMP velocity in the laboratory frame. The 
reduced mass mr,N is defined by mr,N = m^mf^/{m^ + tun), where m-^ is 
the WIMP mass and ton that of the target nucleus. Finally, Umin = a^/Q 
with a = y^TON/STO^j^ is the minimal incoming velocity of incident WIMPs 
that can deposit the energy Q in the detector, and Vmax is related to the 
escape velocity from our Galaxy at the position of the Solar system. 

The total WIMP-nucleus cross section cto in Eq.(l) depends on the 
nature of the WIMP couplings on nucleons. Generally speaking, one has 
to distinguish spin-independent (SI) and spin-dependent (SD) couplings. 
Through e.g., squark and Higgs exchanges with quarks, Majorana WIMPs 
e.g., neutralinos in the supersymmetric models, can have a SI scalar inter- 
action with nuclei ' : 



^Si 



r,p 



7^' 



7^' 
XP 



- ] <pi/pr 



(2) 



where A is the atomic number of target nucleus, TO^.p is the reduced mass 
of WIMPs and protons, and /p is the effective XXPP four-point coupling. 
Note here that the approximation /„ ~ /p predicted in most theoretical 
models has been adopted and the tiny mass difference between a proton 
and a neutron has been neglected. 

Meanwhile, through e.g., squark and Z boson exchanges with quarks, 
WIMPs can couple to the spin of the target nuclei. The total cross section 

1 9 

for the spin coupling can be expressed as ' 



^0 



32 



'-r,N 



J+1 



(5'p)ap + (S'n)an 



(3a) 



and 



^X(p,n) 



^(P,n)l 



(3b) 



Here Gp is the Fermi constant, J is the total spin of the target nucleus, 

(5'(p,n)) are the expectation values of the proton and the neutron group 
spins, and a(p_n) is the effective SD WIMP coupling to protons and neutrons. 
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2. Determining the WIMP mass 



It has been found that the one dimensional velocity distribution function 
of incident WIMPs, fi{v), can be solved analytically from Eq.(l) directly'' 
and, consequently, its generalized moments can be estimated by^ 

(■i;"')('y(Qniin),w(Qmax)) 
/■f(Qmax) 

= / v^h{v)dv 

•'f (Qmin) 



2Q 



("+l)/2, 



'^{Qmin^/F' (^min) H" H" l)-^ri(^min5 ^max) 



max J 



(4) 



Here v{Q) = a\/Q, (5(inin,max) Skve the minimal and maximal cut-off ener- 
gies of the experimental data set, respectively, r{Q^in) = {dR/dQ)Q^Q^^^ 

and 



is an estimated value of the scattering spectrum at Q 
IniQmin, Qmax) Can be estimated through the sum: 



Qr, 



-^n(Qmin) Qmax) — ^ ^ 



(n-l)/2 



(5) 



where the sum runs over all events in the data set between Qmin and Qmax- 
By requiring that the values of a given moment of fi{v) estimated by 
Eq.(4) from two detectors with different target nuclei, X and Y, agree, a 
general expression for determining appearing in the prefactor a" on the 
right-hand side of Eq.(4) has been found as^: 



m 



y/rnxmy 

T^n,x/'R-n,Y - \^mx/mY 



where 
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IQ.X 



(6) 



(7) 



and TZn^Y can be defined analogously. Here n ^ 0, m(^x,Y) and F(^x,y){Q) 
are the masses and the form factors of the nucleus X and Y, respectively. 
Note that, since the general moments of fi{v) estimated by Eq.(4) arc 
independent of the WIMP nucleus cross section ao, the estimator (6) of 
can be used either for SI or for SD scattering. 

Additionally, since in most theoretical models the SI WIMP-nucleus 

cross section given in Eq.(2) dominates^ , and on the right-hand side of 
Eq.(l) is in fact the minus-first moment of fi{v), which can be estimated 
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by Eq.(4) with n = 

PolUf 



-1, one can find that^ 



(: 



TON 



(Qmin) 



lo 



(8) 



Here £ is the exposure of the experiment whieh relates the actual counting 



on 



rate to the normalized rate in Eq.(l). Since the unknown factor pol/p 
the left-hand side above is identical for different targets, it leads to a second 
expression for determining m^^: 



'^a,x/T^a,Y - (mx/mY)^^'^ 
Here m(x,F) oc has been assumed, 



T^a,x = -pr- 

c-x 



'^Qmhi,x'>^x{Qmin,x) 



+ Io,x 



min.X ) 



(9) 



(10) 



and similarly for TZa,Y- 

In order to yield the best-fit WIMP mass as well as its statistical error 
by combining the estimators for different n in Eq.(6) with each other and 
with the estimator in Eq.(9), a function has been introduced*^ 



(11) 



where 



fi,X = 
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axTZi 



X 



300 km/s 



fun 



.+1,X 



for i 



^x 



a,X 



-1, 1, 2, 



mx 



(12a) 



(12b) 



the other rimax + 2 functions /i,Y can be defined analogously. Here rimax 
determines the highest moment of fi{v) that is included in the fit. The fi 
are normalized such that they are dimensionless and very roughly of order 
unity. Note that the first rimax + 1 fit functions depend on through the 
overall factor a and that in Eqs.(12a) and (12b) is now a fit parameter, 
which may differ from the true value of the WIMP mass. Moreover, C is 
the total covariance matrix. Since the X and Y quantities are statistically 
completely independent, C can be written as a sum of two terms: 



Cij = cov {fi,x,fj,x) + cov (/i,r, fj,Y) 



(13) 
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Finally, since we require that, from two experiments with different target 
nuclei, the values of a given moment of the WIMP velocity distribution 
estimated by Eq.(4) should agree, this means that the upper cuts on fi{v) 
in two data sets should be (approximately) equal^ This requires that^ 

Qmax.X • (14) 

ay J 

Note that a is a function of the true WIMP mass. Thus this relation for 
matching optimal cut-off energies can be used only if is already known. 
One possibility to overcome this problem is to fix the cut-off energy of the 
experiment with the heavier target, minimize the x^if^x) function defined 
in Eq.(ll), and estimate the cut-off energy for the lighter nucleus by Eq.(14) 
algorithmically" . 

As demonstration we show some numerical results for the reconstructed 
WIMP mass based on Monte Carlo simulations. The upper and lower 
bounds on the reconstructed WIMP mass are estimated from the require- 
ment that exceeds its minimum by 1. ^*Si and ^^Ge have been chosen as 
two target nuclei. The scattering cross section has been assumed to be dom- 
inated by spin-independent interactions. The theoretically predicted recoil 
spectrum for the shifted Maxwellian velocity distribution (vq = 220 km/s, 
Ve = 231 km/s) ' ' ^ with the Woods-Saxon elastic form factor"' ^' have 
been used. The threshold energies of two experiments have been assumed to 
be negligible and the maximal cut-off energies are set as 100 kcV. 2 x 5,000 
experiments with 50 events on average before cuts from each experiment 
have been simulated. In order to avoid large contributions from very few 
events in the high energy range to the higher moments^, only the moments 
up to rijiiax = 2 have been included in the fit. 

In Fig. 1 the dotted (green) curves show the median reconstructed 
WIMP mass and its la upper and lower bounds for the case that both 
Qmax have been fixed to 100 kcV. This causes a systematic underestimate 
of the reconstructed WIMP mass for input WIMP masses ^ 100 GeV'''. 
The solid (black) curves have been obtained by using Eq.(14) for match- 
ing the cut-off energy Qmax,Si perfectly with Qmax.Go = 100 keV and the 
true (input) WIMP mass, whereas the dashed (red) curves show the case 
that Qmax.Ge = 100 keV, and (5max,Si has been determined by minimiz- 
ing x^(mj^). As shown here, with only 50 events from one experiment, the 
algorithmic process seems already to work pretty well for WIMP masses 



tHere the threshold energies have been assumed to be neghgibly small. 
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1000 



Fig. 1. Results for the median reconstructed WIMP mass as well as its Itr statistical 
error interval based on the x^^fit in Eq.(ll). See the text for further details. 



up to ~ 500 GeV. Though for ^100 GeV determined in tliis way 
overestimates its true value by 15 to 20%, the true WIMP mass always 
lies within the median limits of the la statistical error interval up to even 



= 1 TeV 



6 



On the other hand, in order to study the statistical fluctuation of the 
reconstructed WIMP mass by algorithmic Qmax matching in the simulated 
experiments, an estimator 5m has been introduced as^ 



dm ■ 



1 



if 



< m 



X.lol 



if ra^,\oi < my^,in < rn^,rec ; 
if 



(15) 



if 



Here m^jn is the true (input) WIMP mass, m^^rcc its recon- 
structed value, m--^ ioi(2) are the 1 (2) a lower bounds satisfying 
X^(nT-x,io(i,2)) = X^("^x,rcc) + 1 (4), and m^,hii(2) are the corresponding 
1 (2) a upper bounds. 

Figures 2 show the distribution of the estimator Sm calculated from 
5,000 simulated experiments for WIMP masses of 50 GeV (left) and 200 
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Fig. 2. Normalized distribution of the estimator 5m defined in Eq. (15) for WIMP 
masses of 50 GeV (left) and 200 GeV (right). Parameters and notations are as in Fig. 1. 
Here the bins at Sm = ±5 are overfiow bins, i.e., they also contain all experiments with 
\5m\ > 5. 



GeV (right). For the hghter WIMP mass, simply fixing both Qmax values 
to 100 keV still works fine (the dotted (green) curves in Fig. 1). However, 
the distributions for both fixed Qmax and optimal Qmax matching show 
already an asymmetry of the statistical uncertainties with m^^hii — w^^rcc > 
W;^,rcc — rn^^ioi- The overestimate of light WIMP masses reconstructed by 
algorithmic Qmax matching shown in Fig. 1 is also reflected by the dashed 
(red) histogram here, which has significantly more entries at positive values 
than at negative values. Moreover, these distributions also indicate that 
the statistical uncertainties estimated by minimizing x^(mj^,) are in fact 
overestimated, since nearly 90% of the simulated experiments have \Sm\ < 
1^, much more than ^ 68% of the experiments, which a usual Icr error 
interval should contain. 

For the heavier WIMP mass of 200 GeV, as shown in the right frame 
of Figs. 2, the situation becomes less favorable. While the distributions for 
both fixed Qmax and optimal Qmax matching look more non-Gaussian but 
more concentrated on the median values, the distribution for algorithmic 
Qmax matching spreads out in the range — 1 < Sm < 2. It has even been 
observed that, for larger samples (e.g., with 500 events on average) the 
outspread distribution becomes broader^. Hence, the statistical fluctuation 
by the algorithmic procedure for determining Qmax of the experiment with 
the lighter target nucleus by minimizing could be problematic for the 
determination of if WIMPs are heavy. 
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3. Estimating the SI WIMP— proton coupling 

As shown in the previous section, by combining two experimental data sets, 
one can estimate the WIMP mass m^^ without knowing the WIMP-nucleus 
cross section ctq- Conversely, by using Eq.(8), one can also estimate the 
SI WIMP-proton coupling, |/pp, from experimental data directly without 
knowing the WIMP mass^. 

In Eq.(8) the WIMP mass on the right-hand side can be determined 
by the method described in Sec. 2, r{Qniin) and Iq can also be estimated 
from one of the two data sets used for determining or from a third ex- 
periment. However, due to the degeneracy between the local WIMP density 
Po and the coupling |/pp, one cannot estimate each one of them without 
making some assumptions. The simplest way is making an assumption for 
the local WIMP density po. 

Figure 3 shows the reconstructed SI WIMP-proton coupling as a func- 
tion of the input WIMP mass.. The WIMP mass has again been recon- 
structed with ^*Si and ''^Ge. In order to avoid complicated calculations of 
the correlation between the error on the reconstructed and that on the 
estimator of /q, a second, independent data set with Ge has been chosen as 
the third target for estimating Iq. Parameters are as in Fig. 1, except that 
the SI WIMP-proton cross section has been set as 10~^ pb. 



o 



I CE^ctn un rm 1 



10 50 100 500 

mchi in[GeVl 

Fig. 3. The reconstructed SI WIMP-proton coupling as a function of the input WIMP 
mass. The (red) squares indicate the input WIMP masses and the true values of the 
coupling. The (blue) circles and the (blue) crosses indicate the reconstructed couplings 
and their Irr statistical errors. Parameters are as in Fig. 1 
as 10~* pb. See the text for further details. 



in addition (t|p has been set 
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It can be seen that the reconstructed |/pp arc underestimated for WIMP 
masses ^ 100 GeV. This systematic deviation is caused mainly by the 
underestimate of 7o- However, in spite of this systematic deviation the true 
value of l/pp always lies within the la statistical error interval. Moreover, 
for a WIMP mass of 100 GcV, one could in principle already estimate the SI 
WIMP-proton coupling with a statistical uncertainty of only ~ 15% with 
just 50 events from each experiment. Note that this is much smaller than 
the systematic uncertainty of the local Dark Matter density (of a factor of 
2 or even larger). 

4. Determining ratios of WIMP— nucleon cross sections 

So far we have discussed only the case that the spin-independent WIMP- 

nuclcus interaction dominates. In this section we turn to consider the case 
of the spin-dependent cross section as well as of a general combination of 
these two cross sections. 

4.1. Determining the a^/a^ ratio 

Consider at first the case that the SD WIMP -nucleus interaction dominates. 
By substituting tTg^ in Eq.(3a) and {v~^) estimated by Eq.(4) into Eq.(l) 
and combining two data sets with different target nuclei, an expression for 
the ratio between two SD WIMP-nucleon couplings can be given as 

/an\^° _ {Sp)xT^J,n,Y ± {Sp)Y'R'J,n,X ^-^^^ 
\^pJ±,n {Su)xT^J,n,Y ± ('S'n)y'7?.J,n,X ' 

with 

1/2 

(17) 

and similarly for TZj^n,Y, where n 7^ 0. Note that an/ftp can be estimated 
from experimental data directly through estimating TZn,Xj T^a,x and two 
Y terms by Eqs.(7) and (lO)* without knowing the WIMP mass. 

Because the couplings in Eq.(3a) arc squared, we have two solutions 
for a-a/ap here; if exact "theory" values for 'R'j^n,(x,Y) arc taken, these 
solutions coincide for an/flp = —{Sp)x/{Sn)x and —{Sp)Y/{Sn)Y, which 
depends only on the properties of target nuclei^. Moreover, one of these 



tNotc that the form factor F^{Q) here must be chosen for the SD cross section. 

S Some relevant spin values of the nuclei used for our simulations shown in this paper are 

given in Table. 1. 



n. 



J,n,X 



Jx 



Tin 



Jx + l TlnJ 
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Table 1. List of the relevant spin values of the nuclei used for simulations shown in 
this paper (Data from Ref. 10). 



nucleus 


Z 


J 






-(5p)/(S„> 




8 


5/2 





0.495 





23Na 


11 


3/2 


0.248 


0.020 


-12.40 


37C1 


17 


3/2 


-0.058 


0.050 


1.16 


73 Ge 


32 


9/2 


0.030 


0.378 


-0.079 



two solutions has a pole at the middle of two intersections, which depends 
simply on the signs of {Sn)x and {Sn)Y- since TZj^n.x and TZj.n.Y are always 
positive, if both of {Sn)x and (S'n)Y are positive or negative, the "minus" 
solution (an/ap)^°„ will diverge and the "plus" solution (an/ap)+^„ will 
be the "inside" solution, which has a smaller statistical uncertainty (see 
Figs. 4); in contrast, if the signs of {Sn)x and (S'n)y are opposite, the 
"minus" solution (an/ap)^^„ will be the "inside" solution. 

Figures 4 show the reconstructed (an/flp)^^ estimated by Eq.(16) with 
n = 1 as functions of the true (input) a^/a^ for a WIMP mass of 100 GeV 
(left) and as functions of the input WIMP masses for a^/a-p = 0.7 (right), 
respectively. The shifted Maxwellian velocity distribution with a form factor 
calculated in the thin-shell approximation for the SD cross section^' 
has been used. Parameters are as earlier, except that the minimal cut-off 
energy has been increased to 5 keV for both experiments. Here we have 



Qmm ?5 keV, Qma» ^ 100 ke 



hi = 100 GeV, Ge73 + CI37, 2 » 50 e> 



S8* 



0.5 



■ an/ap[rGG, 1] 
- an/ap[rec, 1, sli] 



Qmm > 5 keV. Qmai: <. 100 keV. an/ap = 7, Ge73 + CI37,2 v 50 
an/ap[reG, 1] 

; =„;sp[rt., 1, sh] 



10 20 30 50 70 100 200 300 300 7001000 

mchilinl [QeVl 



Fig. 4. Preliminary results for the reconstructed (an/cip)^'^ estimated by Eq.(16) with 
n = 1 as functions of the true (input) an/flp (left frame, for a WIMP mass of 100 
GeV) and as functions of the input WIMP mass rrt^ (right frame, for a^/ap = 0.7), 
respectively. See the text for further details. 
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chosen ^^Gc and ^^Cl as two target nnclei in order to test the range of 
interest < a^/a-p < l^"^' and (an/ap)^y_ ^ are thus the inside/outside 
solutions. 

For estimating the statistical error on Un/ap, one needs to estimate the 
counting rate at the threshold energy, r((3min), and its statistical error, 
C7(r((5min))- It has been found that, instead of Qminj one can estimate the 
counting rate and its statistical error at the shifted point Qs,i (from the 
central point of the first bin, Qi)^: 



Qs,i — Qi 



ki 



In 



sinh(fcifei/2) 



kibi/2 



(18) 



where ki is the logarithmic slope of the reconstructed recoil spectrum in the 
first Q-bin and bi is the bin width. We see in the right frame of Figs. 4 very 
clearly that, for WIMP masses ^30 GeV, the la statistical error estimated 
with Qs,i (the dash-dotted (red) lines, labeled with "sh") is ~ 7%, only 
1/3 of the error estimated with Qmin (the dashed (blue) lines). 

One more advantage with using Qs,i instead of Qmin is that the sta- 
tistical error on a^/cip estimated with different n (namely with different 
moments of the WIMP velocity distribution) at Q = Qs,i are almost equal. 
Therefore, since 



J-l,X = 



Jx 



2rxiQx,sA) 



Jx + lJ SxFliQx.ss] 



1/2 



(19) 



one needs thus only events in the low energy range (~ 20 events between 5 
and 15 keV in our simulations) for estimating Un/ap. 



4.2. Determining the '''^p/n/'^xp 

Now let us combine WIMP nucleus scattering induced by both SI and 
SD interactions given in Eqs.(2) and (3a) (with the corresponding form 
factors). By modifying F^{Q) and /„ in the estimator (4) of the moments 
of the WIMP velocity distribution, the ratio of the SD WIMP-proton cross 
section to the SI one can be solved analytically as^' I' 

^XP Fgj^^{Q^in,x)T^m,Y — FgjY{Qmin,Y)'^m,X ^^q^ 

C^p Cp^xF^-Q-^{Qmin,x)Tlm,,Y — Cp,r ^SD,y (Omin,l')7^m,X 



^In this section we consider only the case of f^p , but all formulae given here can be 

applied straightforwardly to the case of by exchanging n ^ p. 
"Qmin appearing in this section can be replaced by Qs,i everywhere. 
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Here 



and 



_ 4 ( Jx + l 



{Sp)x + {Sn)x{an/ap 



iX 



(21) 



(22) 



TZm,Y and Cp^y can be defined analogously. Note that a "minus (— )" sign 
appears in tlie expression (20). 

By introducing a third target having only the SI interaction with 
WIMPs, ttn/ap appearing in Cp^x and Cp^y can again be solved analyti- 
cally as 

SI-|-SD / \ /^^^^—^^^^— I 

_ - {Cp,XSn/p,X - Cp,ySn/p,rj ± ^/Cp^xCp^y Sn/p,X " •Sn/p,F | 



Or 



Here 



Cp,X = ^ 



2 2 
^P>-^'^n/p,X ~ ^P'YSn/p,Y 



_ 4 ^Jx + l^ f{Sp)x\ ^2 \ 



(23) 



S\ Jx 



;i2 N / ^m,r \ j-,2 



-^SI.zCQmin.z) 



7e. 



~ -Psi,v(<3min,F) 



(24) 



Cp,!' can be obtained by simply exchanging X <-> F, and Sn/p = {Sn)/{Sp). 
Howcivcr, in order to reduce the statistical uncertainties contributed from 
estimate of On/flp involved in Cp^x and Cp^y, one can use one target with 
the SD sensitivity (almost) only to protons or to neutrons combined with 
another one with only the SI sensitivity. For this case Cp^x is independent 
of Un/ap and the expression (20) for a^/a': 



can be reduced to 



^SD 

XP 

XP 



n 



■Y ~ Fgj y{Qthre,Y)^rn,X 



-'P.^-^SD,X 



m.Y 



(25) 



Figures 5 show the reconstructed (^^p/c^jE) (left) and a^/a^p (right) 
estimated by Eqs.(20) and (25) as functions of the true (input) a^/ap, 
respectively. Besides ''^Ge and ^''Cl, ^^Si has been chosen as the third target 
for estimating Un/op by Eq.(23); whereas ^^Ge has been chosen as the 
second target having only the SI interaction with WIMPs and combined 
with 23Na (for cr^p/^xp) (^^^ ^^I'^xp) ^^ing Eq.(25). We see 

here that, since the SD WIMP-nucleus interaction doesn't dominate for 
our simulation setup, f^p /c^p estimated by Eq.(2G) has two discontinuities 
around the intersections at a^/ap = —0.079 and especially at On/flp = 1.16, 
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sigmaSIp = 1D"(-B) pb, ap = 1 , mchi = 100 GeV, 2!?j « 50 events 5 - 100 keV, sigmaSIp = 10"(-B) pb, ap = 1 , mchi = 100 GeV, 2rd « 50 ei 
■ Gc73 + CB7+ Si28 



■ Ge76 + Na23 



Gc73 + C137 + Si28 
■ Gc76 + 017 



/ 

// 



--^^ . 1 J //JO' 



-1 -0.5 0.3 1 1.5 2 -1 -0.5 " 0.5 1 1.5 2 

an/apfin] an/apfin] 

Fig. 5. Preliminary results for the reconstructed cr^p / tr^p (left) and tr^J^ / cr^p (right) as 
functions of the true (input) an /flp , respectively. The dashed (blue) curves indicate the 
values estimated by Eq.(20) with an/flp estimated by Eq.(23); whereas the dash— dotted 
(red) curves indicate the values estimated by Eqs.(25). cr^p and a-p have been set as 10~* 
pb and 0.1, respectively. The other parameters are as in Figs. 4. Note that, since we fix 
(T^p and Op, (T^p/crlp shown here is a constant, whereas cr^J^/cr^p oc a parabola. 



the intersection determined by the —{Sp) / (Sn) value of '^^Cl. However, from 
two experiments with only ~ 20 events in the low energy range, one could 
in principle already estimate f^p/f^p and a^/a^p by using Eq.(25) with 
a statistical uncertainties of ^ 35%. 



5. Summary and conclusions 

In this article we described model-independent methods for determining the 
WIMP mass and their couplings on nucleons by using future experimental 
data from direct Dark Matter detection. The main focus is how well we 
could extract the nature of WIMPs with positive signals and which problems 
we could meet by applying these methods to (real) data analysis. 

In Sees. 2 and 3 we discussed the determinations of the WIMP mass and 
its spin-independent coupling on protons. If WIMPs are light (m^ ~ 50 
GeV), with 0(50) events from one experiment, their mass and SI coupling 
could be estimated with errors of ~ 35% and ^ 15%, respectively. However, 
in case WIMPs are heavy (m,^ ^ 200 GeV), the statistical fluctuation by 
the algorithmic procedure for matching the maximal cut-off energies of the 
experiments could be problematic for estimating their mass, and thereby 
their SI coupling. 

In Sec. 4 we turned to consider the spin-dependent interaction. The 
simulations show pretty small statistical uncertainties. Moreover, differ- 
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ing from the traditional method for constraining the SD WIMP-nucleon 
couplings^^' we do not make any assumptions on po, fi{v), and 

ru;^. The price one has to pay for this is that positive signals in at least 
two different data sets with different target nuclei are required. In addi- 
tion, without independent knowledge of po, one can only determine ratios 
of cross sections. 

In summary, once two (or more) experiments measure WIMP events, 
the methods presented here could in principle help us to extract the nature 
of halo WIMPs. This information will allow us not only to constrain the 
parameter space in different extensions of the Standard Model, but also to 
confirm or exclude some candidates for WIMP Dark Matter 

Acknowledgments 

This work was partially supported by the Marie Curie Training Research 
Network "UniverseNet" under contract no. MRTN-CT-2006-035863, by 
the European Network of Theoretical Astroparticle Physics ENTApP IL- 
IAS/N6 under contract no. RII3-CT-2004-506222, as weU as by the BK21 
Frontier Physics Research Division under project no. BA06A1102 of Korea 
Research Foundation. 

References 

1. G. Jungman, M. Kamionkowski, and K. Griest, Phys. Rep. 267, 195 (1996). 

2. G. Bertone, D. Hooper, and J. Silk, Phys. Rep. 405, 279 (2005). 

3. P. F. Smith and J. D. Lewin, Phys. Rep. 187, 203 (1990). 

4. J. D. Lewin and P. F. Smith, Astropart. Phys. 6, 87 (1996). 

5. M. Drccs and C. L. Shan, J. Cosmol. Astropart. Phys. 0706, Oil (2007). 

6. M. Drees and C. L. Shan, J. Cosmol. Astropart. Phys. 0806, 012 (2008). 

7. C. L. Shan and M. Drees, arXiv: 0710. 4296 [hep-ph] (2007). 

8. J. Engel, Phys. Lett. B 264, 114 (1991). 

9. M. Drees and C. L. Shan, arXiv: 0809. 2441 [hep-ph] (2008). 

10. F. Giuliani and T. A. Girard, Phys. Rev. D 71, 123503 (2005). 

11. H. V. Klapdor-Kleingrothaus, I. V. Krivosheina, and C. Tomei, Phys. Lett. 
B 609, 226 (2005). 

12. V. A. Bcdnyakov, Phys. Atom. Nucl. 67, 1931 (2004). 

13. J. Ellis, K. A. Olive, and C. Savage, Phys. Rev. D 77, 065026 (2008). 

14. D. R. Tovey et ai, Phys. Lett. B 488, 17 (2000). 

15. F. Giuliani, Phys. Rev. Lett. 93, 161301 (2004). 

16. T. A. Girard and F. Giuliani, Phys. Rev. D 75, 043512 (2007). 

17. G. Bertone et al, Phys. Rev. Lett. 99, 151301 (2007). 

18. V. Barger, W. Y. Keung, and G. Shaughnessy, Phys. Rev. D 78, 056007 
(2008). 



